New interactions in the dark sector mediated by dark energy 
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Cosmological observations have revealed the existence of a dark matter sector, which is commonly 
assumed to be made up of one particle species only. However, this sector might be more complicated 
than we currently believe: there might be more than one dark matter species (for example two 
components of cold dark matter or a mixture of hot and cold dark matter) and there may be new 
interactions between these particles. In this paper we study the possibility of multiple dark matter 
species and interactions mediated by a dark energy field. We study both the background and the 
perturbation evolution in these scenarios. We find that the background evolution of a system of 
multiple dark matter particles (with constant couplings) mimics a single fluid with a time-varying 
coupling parameter. However, this is no longer true on the perturbative level. We study the case of 
attractive and repulsive forces as well as a mixture of cold and hot dark matter particles. 



I. INTRODUCTION 

Cosmological observations of the largest structures in 
the universe, the anisotropies in the cosmic microwave 
background (CMB) radiation and Type la supernovae 
suggest that normal matter contributes only 4% of the 
total energy budget of the universe [1, 2, 3]. The major- 
ity of energy and matter is made up of cold dark matter 
(CDM, 21%), which clusters on scales of galaxies and 
clusters of galaxies, and a much more homogeneous dis- 
tributed dark energy component (75%). Whereas dark 
matter candidates are well motivated within particle the- 
ories beyond the standard model, it is fair to say that 
the discovery of the accelerated expansion of the uni- 
verse came as a surprise and that it will be much harder 
to interpret within particle physics theories. Candidates 
for dark energy include the cosmological constant (e.g. 
[4]), scalar fields (see e.g. [5, 6, 7, 8, 9, 10, 11, 12, 13, 
14, 15, 16, 17, 18, 19, 20, 21]) or vector fields (see e.g. 
[22, 23, 24]). When analysing the data it is usually as- 
sumed that General Relativity is the correct theory of 
gravity at the largest visible scales. Of course, dark en- 
ergy might signal a breakdown of General Relativity it- 
self, a possibility which has been recently studied in detail 
(see e.g. [25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 
38, 39, 40, 41] and references therein for recent work). 
In this paper, however, we will assume that dark energy 
can be attributed to a slow-rolling scalar field, interact- 
ing with other matter forms present in the universe, most 
notably with CDM and hot dark matter (HDM, in the 
form of massive neutrinos for example). 

Since the nature of dark energy remains unknown, it is 
important to study new effects resulting from dark energy 
and its possible interactions with the rest of the world. 



*Email address: php04awb@shcfBcld.ac.uk 
tEmail address: C.vandcBruck@shcfficld.ac.uk 
tEmail address: Lisa.Hall@shcfHcld.ac.uk 



Whereas new forces between normal matter particles are 
heavily constrained by observations (e.g. in the solar sys- 
tem and gravitational experiments on Earth), this is not 
the case for neutrinos and dark matter particles. At the 
moment, constraints on new forces between neutrinos or 
dark matter can only come from cosmological observa- 
tions. Additionally, it is usually assumed that there is 
only one type of dark matter; this is well motivated by 
supersymmetrical models which predict that dark matter 
is the lightest supersymmetric particle. However, there 
might be more than one dark matter species. Such a 
possibility has been investigated previously in [19] and 
can be motivated from string theory [42, 43]. In this 
paper, we study the cosmological implications of a dark 
sector with multiple matter species (CDM, HDM) and 
consider the possibility that new forces are mediated be- 
tween the particles by a dark energy scalar field. The 
aim is to study both the background evolution as well as 
the evolution of perturbations in the CDM and/or HDM 
fluids. We take into account the full evolution of the 
HDM component, both for the background as well as on 
the perturbative level [44, 45]. We will see that there 
are new effects associated with new forces in the CDM 
and/or HDM sector, and there is much hope that cos- 
mological data will strongly constrain such interactions. 
We do not address the question of whether the effects of 
the new interactions arc inherent in alternative theories 
of gravity, a possibility which has been pointed out be- 
fore [46] . In both coupled quintessence (including models 
considered here) and modified gravity, the growth rate of 
density fluctuations will be modified due to the presence 
of new degrees of freedom. This may lead to a degener- 
acy between these models when only linear perturbations 
are taken into account. 

The new force mediated by the dark energy scalar field 
is always attractive for particles of the same kind. How- 
ever, if the couplings for the individual species (to be de- 
fined in the next section) have opposite signs, the force 
between different forms of matter will be repulsive. For 
the cases studied in this paper the force mediated by dark 
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energy is always long-ranged. As a result, the total force 
acting on the dark sector particles (i.e. the sum of the 
gravitational force and the new force) can be smaller than 
the gravitational force itself or even repulsive, depending 
on the coupling. In addition, we will see that when the 
couplings have opposite signs, the expansion of the uni- 
verse in the matter dominated era can behave as if there 
is no new interaction between the particles. Cosmological 
structure formation, however, is affected by the presence 
of the new forces. 

When studying concrete examples, we will assume that 
there are only two dark species and one dark energy 
scalar field. The paper is organized as follows. In the next 
Section, we describe the background evolution for three 
typical models. In Section 3, we study the perturbations 
in these models, calculating both CMB anisotropy and 
matter power spectra. We conclude in Section 4. 



II. BACKGROUND EVOLUTION 

We consider a system of N species of matter coupled 
to a single scalar field </>, which we assume will be re- 
sponsible for the late-time accelerated expansion of the 
universe. We assume a flat, homogeneous Friedmann- 
Robertson- Walker universe with line-element given by 



ds 2 = a{r) 2 (-dr 2 + 



S t jdx l dx J ) 



(1) 



The Einstein equations, describing the evolution of the 
scale factor a(r), are given by 



H 2 = 
H = 



a 
a 

4irGa 2 



P, 



8irGa 2 
(p + 3p) , 



(2) 
(3) 



where p(r) and p(r) are the total energy density and 
pressure respectively, and overdots represent differentia- 
tion with respect to conformal time r. In the following 
we set 87rG = 1 . 

As is customary, we define /3j as the strength of the 
coupling between the i th matter species and the scalar 
field, which is given by 



0i 



din rrii((f>) 



(4) 



where m, is the particle mass of species i. Conservation of 
energy momentum for each of the coupled fluids requires 
that 



which yields 



Pi + 3H(pi + pi) = (3i4>{pi - 3p t ) 



(5) 



(6) 



in the Robertson- Walker spacetime for each of the in- 
dividual coupled matter species. The evolution of the 



scalar field is described by the Klein-Gordon equation, 
which in the presence of matter couplings is given by 



,dV 



N 



2H<p + a 2 — = -a 2 J2&(Pr-3p l )- (7) 



In order to investigate the late-time attractor solutions 
of this system, we perform the usual transformation of 
variables 





u., = — \ — 



with 



N 

2 i 2 , 2 i V"^ 2 

x + y + z + 2ju\ 

i=l 



1. 



(8) 



In these equations, p~ l and pb are the energy density of 
radiation and baryons respectively. Note that in the fol- 
lowing we assume that Pi is constant and Ui is non- 
relativistic (which is valid for late times). We also use 
the Friedmann constraint (Eqn. (8)) to eliminate v and 
consider a general quintessence potential with 



A 



1 dV 
V~d4' 



(9) 



Differentiation with respect to a = In a yields the fol- 
lowing set of equations, 
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(10) 

(ii) 

(12) 
(13) 

(14) 



It has been shown in [16] that for a single coupled fluid 
with an exponential coupling and exponential potential 
there are two late time attractor solutions consistent with 
an accelerating universe. The critical point for the first 
solution (Attractor A) is characterised by fi^ = 1 at late 
times, and leads to acceleration provided A < s/2. The 
second critical point (Attractor B) is a scaling solution, 
with Qfi dependent upon the choice of f3 and A. This 
solution results in acceleration if A < 2(3. For full details 
we refer to [16]. 
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Case 


A 


Species 


Coupling 


Relative Densities 


1 


2 


ft 




(&)„ 
i __/_u 


i 


A 


CDM 


CDM 


4.0 


0.1 


5.2 x 10~ 4 


ii 


/3 

V 2 


CDM 


CDM 


-1.5 


l 


0.28 


in 


2 


HDM 


CDM 


5.8 


-0.1 


0.56 



TABLE I: Cases considered in this paper: Two cases (I, II) 
with two species of CDM and one case (III) with CDM and 
HDM. For all cases ^initial = 0.1 and we consider a flat uni- 
verse with cosmological parameters: h = 0.7, fi^ h 2 = 0.022 
and flQQ M /i 2 = 0.12 for the coupled CDM species. 



We find that introducing multiple coupled fluids does 
not affect this result because at late times the system re- 
duces to the single fluid case as the highly coupled fluids 
dilute at a faster rate. Indeed, the only non-trivial ad- 
ditional critical point which arises due to the addition of 
multiple coupled species is characterised by 



N 

E 

»=i 



ftu? = 0, 



(15) 



with x = y = z = and 17^ = 0. This solution does 
not lead to an accelerating cosmology, but can play an 
important role as a transient solution during the matter 
dominated epoch, as we will discuss later. 

In order to obtain an understanding of the background 
evolution, it turns out to be useful to combine the cou- 
pled matter species into a single effective fluid, with total 
energy and pressure given by 



N 

E 



Pi, Pc 



N 

E; 



(16) 



and an effective coupling /3 e e given by 



e 



N 

J2fo {Pi - 3Pi) 

i=l 

Pc ~ 3p c 



(17) 



We can therefore recast our multi-fluid system of conser- 
vation equations into a single conservation equation for 
p c with an effective coupling /3 off . We will see later that 
this is only possible for the background. It should be 
further noted that the effective coupling is a function of 
both the scalar field and the energy densities of the cou- 
pled fluids, and therefore will in general vary in time even 
if the couplings of the individual fluids are constant. In 
other words, an observer, studying the background evo- 
lution and assuming a single fluid only, will deduce that 
the coupling evolves in time, although there are several 
species with constant individual couplings. 

We now consider the effect on the cosmological back- 
ground of coupling multiple fluids to the scalar field. For 
simplicity, we assume that there are two species of cou- 
pled matter, and that the form of the coupling between 



the coupled fluids and the scalar field is of exponential 
form for concreteness: 



mi(<j>) 



(*) 04 d> 



(18) 



We also take a typical quintessence exponential potential 
for the scalar field, namely 



V{4>) = Me~ M , 



(19) 



where M and A are chosen to give late time acceleration 
today. 

It is instructive to consider three cases, each distin- 
guished by either the sign of the coupling or the type of 
dark matter considered. The specifications can be found 
in Table I. In the first two cases (I, II) we consider two 
types of CDM, but while the individual couplings Pi and 
/?2 are both positive in Case I, the couplings have oppo- 
site signs in Case II. For both of these cases we choose 
A and (3i such that the system approaches the acceler- 
ated late-time attractor with 0,^ — 1 (Attractor A). A 
scenario with mass-varying neutrinos (as a single cou- 
pled species) in which the system approaches Attractor 
A has been studied in [21, 44]. As will be apparent later, 
the chosen parameters lead to a consistent background 
cosmology but the predicted matter power spectra differ 
significantly from the ACDM case and are observation- 
ally ruled out. The purpose of this paper, however, is 
to illuminate the physical effects of these models and the 
extreme parameter choices reflect this. 

It is possible to consider a scenario in which a sys- 
tem with two CDM species reaches the late-time accel- 
erated scaling solution (Attractor B). In order to achieve 
an observationally consistent cosmology, this would re- 
quire one species to remain sub-dominant until today, 
when it would begin to scale with dark energy. In this 
case, the mass of such a species would vary greatly, most 
likely leading it to become relativistic in the past. This is 
the third scenario we consider (Case III), in which a dom- 
inant CDM component and a subdominant HDM com- 
ponent each couple to dark energy 

A similar system, albeit with the CDM uncoupled to 
dark energy, has recently been considered in [47] in the 
context of neutrinos with a growing mass. We choose A 
smaller than in this latter paper. Additionally, we choose 
to take couplings with opposite signs, such that the mass 
of the HDM particle grows in time. Note that for the 
HDM species, we fully consider both the non-relativistic 
and relativistic behaviour of the particles, integrating the 
Fcrmi-Dirac distribution function over phase space [48]. 

In Cases I & II the particles are highly non-relativistic. 
Then, from eqn. (6) and using pi = 0, the evolution of 
the energy density of the i th coupled fluid is given by 



p {°) e /3iO-4>o) 



(20) 



where pf 3 ^ is the energy density of species i at the present 
epoch. With this choice of couplings and two CDM 
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FIG. 1: Relative energy densities for Case 1,11 and III. Solid 
line: Coupled CDM. Long-dashed line: radiation. Short- 
dashed line: dark energy. Dotted line: neutrinos. The faint 
curves show a typical cosmology with /3 = for comparison. 



species, an observer assuming a single fluid would detect 
a time-dependent coupling, which would take the form: 



C 



1 + A 



(21) 



where A = pip\ j p\ , B = p\ — p2, C = P2 and are all 
constant. 

The evolution of the density parameters are shown in 
Fig. 1 for the three cases. The effective coupling experi- 
enced by the combined fluid, as defined in eq. (17), is the 
weighted value of the coupling strengths of the two fluids, 



where the weighting factor is given by (/?, — 3j>i). Hence, 
if the fluid with larger coupling (e.g. species 2 with /3 2 ) 
is dominant at early times, then f5\ < /? c ff < 02 ■ How- 
ever, the species with the higher value of coupling also 
dilutes at a greater rate, and so at later times this species 
becomes sub-dominant, and the effective coupling tends 
towards the lower value for f3. The evolution of {3 e g is 
shown in the left panel of Fig. 2. 

The effect of the couplings on the background evolution 
is to modify the evolution of the energy densities of the 
coupled fluids. Therefore, even if the density parameters 
are tuned to the standard values today (i.e. w 0.7, 
f^CDM ~ 0.25), their values at early times are different 
from what one would expect in the uncoupled case. This 
will affect, among other things, the predictions for the 
CMB anisotropies. 

As can be seen in Fig. 2, for Case II, the effect of 
choosing opposite signs for [3 for the two fluids is to drive 
the effective coupling of the combined fluid towards zero 
during the matter dominated epoch. Therefore for some 
period of time the background cosmology of the system 
behaves essentially identical to the usual uncoupled case. 
We can best understand this behavior by considering a 
transient solution to Eqns (10-14) which is valid when the 
coupled matter species dominates the energy density. In 
this regime we can take u[ = and x 2 ,y 2 , z 2 -c u 2 and 
x',y',z' ~ 0, which leads to a critical point given by 
Eqn (15), which is equivalent to (3 c ff = 0. This transient 
regime ends when the energy density stored in the scalar 
field becomes non-negligible. 

More heuristically, if we instead consider the dynamics 
of the scalar field, then the evolution of <f> is driven by 
the properties of the effective potential 



N 



dV eS dV ^ 



(22) 



which during the matter dominated epoch can be approx- 
imated by 



dV cS 



N 

E 



Pi Pi, 



(23) 



if the potential can be neglected, which is usually the 
case for quintessence potentials and in models where the 
couplings arc of opposite signs. Clearly, in this scenario, 
the effective potential will have a minimum located at 



1 



In 



» 



(24) 



V Pi Pi 

We can expect the field to sit in this minimum if 
4^ > (H/a) 2 , i.e. 



H 



-( Pl +p 2 )<P 2 1 Pi+P 2 2 P2, (25) 
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which requires 

(ft - ft) (ftft + 1/3) > (assuming ft > 0) (26) 

and hence ft > — 3^7- Therefore, provided this con- 
straint is satisfied, and ft and ft are of opposite sign, 
we can expect fta to be driven to zero during the mat- 
ter dominated epoch as the field will be forced to sit in 
the minimum of the effective potential. Note in partic- 
ular that this behavior is independent of the form of V, 
provided that the scalar field satisfies the conditions dis- 
cussed above. When this condition is not satisfied the 
evolution of the system resembles that of Case I. 

Case III is unlike the previous two cases in that the 
values of A and (3 are chosen such that the system ap- 
proaches the accelerated scaling attractor (Attractor B). 
This requirement restricts the choice of density param- 
eter, since tti = fij(A, 0) during the scaling regime (see 
[16]). As the mass of the growing species is expected 
to vary by orders of magnitude we choose to identify 
the mass-growing species with neutrinos, fully consider- 
ing the behavior of the neutrinos in the relativistic and 
non-relativistic regimes. As in Case II, wc choose cou- 
plings of opposite sign, such that ftg is driven towards 
zero during the matter dominated epoch. 

We compare Case III with the growing-mass model 
of [47], who consider the dynamics of Attractor B with 
one coupled species. In our scenario, the field sits in the 
minimum of the effective potential; it becomes heavy due 
to the presence of multiple contributions from the fluids 
with couplings of opposite signs. This is not true in the 
growing-mass model. 

As a final point we consider the evolution of the ap- 
parent equation of state, as discussed in [49]. For a sys- 
tem with multiple species of coupled fluids, the apparent 
equation of state of the dark energy sector as measured 
by an observer assuming an uncoupled system of fluids 
is given by 



where w# = ^ is the equation of state of the scalar field 
and 




Pi 



(0) 



(28) 
(29) 



"'a 



(27) 



It can be seen that in all of these models iu ap is rapidly 
varying and strongly negative, crossing through u> ap = 
— 1 close to 2 ~ 1. Note that in Case III the appar- 
ent equation of state becomes ill-defined, since x goes 
through unity. This is due to our choice of parameters. 



III. PERTURBATIONS 

In the previous section we demonstrated that the back- 
ground evolution of a cosmology with two fluids coupled 
to a single scalar field could be mimicked by a system 
with a single fluid coupled to a scalar field with an ef- 
fective coupling given by eqn. (17). It is therefore im- 
possible to distinguish between the two systems using 
probes of the cosmological background, such as measur- 
ing the expansion history. We have also seen that a sys- 
tem which comprises of multiple fluids with constant cou- 
pling strengths ft may appear to be a single fluid with 
a time varying ft or indeed a system for which j3 = 
for at least some part during the matter dominated era. 
We now investigate the effect that these couplings have 
on the evolution of cosmological perturbations, and de- 
termine whether the degeneracy between the single fluid 
and the multi-fluid approach can be broken. 

For each individual fluid, the equations governing the 
evolution of the density perturbations can be obtained 
from the perturbation of eqn. (5) and, at first order, are 
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given in the synchronous gauge by (see e.g. [50, 51]) 
^ = 3 ( H + prf) ( w 



Sp 



+ Pi (1 - 3wi)5(j) + 



dfa 



(1 - 3wi) 



H(l - 3wi)9i 
1 3Wi -k*56 



k z S, 



1 + Wi 



1 + Wi 1 1 + Wi 

(3i{\ - 3wi)4>9i - k 2 <Ji 



(31) 



where for each species i wc have defined: Si — Spi/pi 
is the density contrast, Spi is the pressure perturbation, 
ft, is a the metric perturbation, 0i is the gradient of the 
velocity field, a is the shear stress and S<p is the pertur- 
bation in the scalar field. For the HDM case, we use the 
full Boltzmann hierarchy to evaluate the density pertur- 
bations [44, 45]. 

The perturbed Klein-Gordon equation, governing the 
evolution of the scalar field perturbation Sep is given 
by [50, 51] 



54> + 2H5(/) + 

N 



2 ' a 2 



d 2 V 



(32) 



dpi 



Pi(5pi - 3Spi) + -T-5(j)(pi - 3pi) 



The evolution of the density contrasts in the different 
cases are shown in Fig. 3, where we consider the limits 
of large and small scales. We now discuss both regimes 
separately. 



Small Scales 



In Case I, where Pi fa > 0, the two species experience 
an additional attractive force between both similar and 
dissimilar particles, increasing the growth of the density 
perturbations. However in Case II, Pi and Pi have oppo- 
site signs, so whilst like particles feel an additional attrac- 
tive force, two dissimilar particles experience a repulsive 
force between each other. This effect can be observed in 
the evolution of the density contrast in Fig. 3(a); Si is 
enhanced in Case I whilst in Case II S2 is enhanced but 
Si is supressed (note that Si becomes negative). For Case 
II, after horizon crossing, the dominant fluid (species 2) 
clusters at an enhanced rate due to the additional contri- 
bution to G cff . Species 1 also feels an additional attrac- 
tive force towards itself, but also experiences a repulsive 
force from species 2. As species 2 dominates the en- 
ergy density, the bulk of the ordinary matter falls into 
the gravitational potential wells that it generates, whilst 
species 1 is repulsed from these overdensitics. 

However, note that during the period when the field 
is in it's minimum, the behavior of the combined fluid 
is identical to that of a single uncoupled species. It is 
only at late times, when the field evolves away from it's 
minimum due to the onset of dark energy domination, 
that the extra coupling force is felt, and the perturbations 
in the total fluid grow at an enhanced rate. 

Case III is similar to Case II, except that the HDM 
undergoes freestreaming until it becomes non-relativistic, 
which suppresses the growth of the density contrast (Si). 



B. Large Scales 

On superhorizon scales the evolution of Si can be ap- 
proximated by 



On small scales, the perturbed Klein Gordon equation 
can be approximated by 



N 



k 2 5(f> + -h4> = -a 2 ^2 fadPi 



(33) 



where we have assumed that pi is constant and Wi ~ 0. 
Using one of the perturbed Einstein equations [48] 



h + Hh = -a 2 ^ Spi 



we can show that 

3 N G off 

St = -(H + pM s\ + -h 2 njS : - 



N 



(34) 



(35) 



where 



Glf = G N (1 + 2pipj) . (36) 
This equation has also been obtained in [51]. 



Oi = — 
2 



■fa 



(37) 



for Wi = 0, which reduces to the usual expression for 
small values of coupling. The suppression of growth 
(prior to horizon crossing, z + 1 > 10 2 ), seen in Fig. 3(b) 
for the highly coupled species, arises due to this addi- 
tional coupling term, as the growth of Si is modified by 
the term PiS(f>. In our models, 8<f> is negative and we ob- 
serve a suppression of 8. However, once the mode enters 
the horizon, Eqn. (37) is no longer valid, the additional 
force between the particles mediated by dark energy can 
be felt and the growth of Si is enhanced. 

We point out that for all of the models considered in 
this paper, the range of the additional scalar force 



x 1 v, 2 d 2 V 
A~_ where m, = ^ 



(38) 



is much greater than the horizon size and we do not there- 
fore observe any effects due to the finite range of the force. 
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z+1 z 1 

(a)Small Scales (b)Large Scales 

FIG. 3: Evolution of the density contrasts for the three cases on small and large scales. Thin solid line: f3 = 0. Short dashed 
line: 5i. Long dashed line: 82- Thick solid line: <5 C . Dotted line (Cases I and II only): equivalent single fluid with j3 — /3 e g- 



C. Comparison to Equivalent Single Fluid 

As for the background, we can consider the two species 
of matter to be a single fluid with 

N ?>P 
Sp c = S pi, 6 C = — -, (39) 

N 

{p c +Pc)dc = ^2{p l +Pi)0 l . (40) 
1=1 



As explained in Section II, an observer studying the back- 
ground evolution, who assumes only a single fluid dark 
matter sector, will infer that the coupling evolves in time, 
Asff = fleff (</>)■ With two dark matter species, an observer 
might detect a coupling of the form given in Eqn. (21). 

However, it is apparent that we cannot re- write Eqns 
(30) and (32) into an equivalent single fluid form, due 
to the ^ term. Whilst for the single fluids (3 is purely a 
function of 0, in the combined fluid (3 becomes a function 
of 4> and pi. The term 8(3 receives contributions from the 
variation of <j) and the variation of Spi, which itself is a 
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FIG. 4: CMB anisotropy spectra and matter power spectra for the three cases. Solid line: (3 = 0. Short dashed line: 2 coupled 
fluids. Long dashed line (Cases I and II only): single fluid with j3 = f3 e g . 



function of <fi and h. Hence, using /3 e(4>) (inferred from 
background observations) to calculate the perturbation 
evolution would lead to incorrect predictions for 6. For 
example, assuming a single fluid with an effective cou- 
pling leads to a different evolution of perturbations than 
if the observer had correctly identified two species. 

In order to highlight this observational difference (al- 
lowing us to break the degeneracy), we calculate the 
evolution of perturbations for an equivalent single fluid 
model with field- varying coupling /3 e ff (0) ■ This is plot- 
ted in Fig. 3 for a direct comparison for Cases I and II. 



(This effect is most apparent when we consider the mat- 
ter power spectrum in the next section.) We also plot 
the case (3 = for comparison. 

In Case II, when (3 e g = 0, the density contrast of the 
combined fluid, 6 C , mimics a single uncoupled fluid, which 
of course coincides with the equivalent single fluid (with 
/3eff = 0) during this time. 

For Case III, however, we do not compare the two fluid 
case with an equivalent single fluid. This is because /3 e ff 
is a function of both the pressure and energy density of 
the coupled species, one of which is highly relativistic 
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at early times. Although (3 e s could be reconstructed as 
/3 e jf (</>), it would not take the simple analytic form given 
in Eqn. (21). 



D. Observational Signatures 

The differing behaviour of the perturbations for the 
multiple fluid and the single effective fluid leaves a de- 
tectable imprint on both the anisotropy spectrum and the 
matter power spectrum, as shown in Fig. 4. Once again, 
we compare the results for the multiple coupled fluids 
with the equivalent single coupled fluid with /3 = p c s(4>). 
We remind the reader that the parameters are chosen 
such that the new effects of the interactions are empha- 
sized. 

For both Cases I and II, the matter power spectra are 
relatively unchanged on large scales, as the growth of 
density perturbations is weakly affected by the coupling 
(see Fig. 3(a)). On small scales, however, the enhanced 
growth of S c is seen explicitly in the matter power spec- 
tra. The difference between multiple and single field flu- 
ids is apparent and breaks the degeneracy seen in the 
background solution. 

In Case III, we see a suppression of power in the mat- 
ter spectrum compared with the standard case. This is a 
consequence of the mass variation of the neutrinos, which 
are lighter in the past than in the standard case. They 
therefore undergo frcestrcaming for longer, which sup- 
press the growth of density perturbations. This effect is 
partly compensated by the increased effective Newton's 
constant, , experienced by the CDM particles, which 
would otherwise clump more than in the standard case. 

Fig. 4(b) shows the effect of the coupled fluids on the 
temperature anisotropy spectrum. It has been shown 
previously that for a single coupled species of CDM, the 
coupling is constrained to \j3\ < 0.1 [52], due to the 
large contribution from the coupling to the late-time In- 
tegrated Sachs-Wolfe (ISW) effect. For Cases I and II, 
however, we observe a reduction in the ISW effect, even 
though we have coupling strengths 0(1). This is because 
the effective coupling strength of the combined fluid in 
both Case I and Case II has been chosen to be 0(0. 1) 
at late times, to ensure that the ISW effect docs not 
contribute too greatly on large scales. It should also be 
noted that, in a system with multiple coupled fluids, the 
fluid with the smallest (3 dominates the energy density 
of the coupled species at late times, as this species di- 
lutes slower than species with larger couplings. In this 
scenario, therefore, we are not restictcd by previous con- 
straints on (3. 

For Case III, we can see that the contribution to the 
ISW from the coupling is enhanced relative to the height 
of the first acoustic peak. 

An additional signature can be found in comparing the 
relative heights of the second and third peaks, which are 
significantly altered in Case I. This arises partly due to 
the change in relative energy densities at the time of last 



scattering as seen in Fig. 1. For Cases II and III, the 
relative energy densities at this time remain unchanged 
as /3 e ff — ► and the fluids mimic standard cosmology. 

We can also see that there is a significant variation in 
the anisotropy spectrum and the matter power spectrum 
for the two coupled species of matter in comparison with 
the equivalent single species of matter. This difference 
therefore leaves a detectable signature, which allows dif- 
ferentiation between coupled fluids of single or multiple 
species. 



IV. CONCLUSIONS 

In this paper, we have considered the effect of new 
forces mediated by dark energy, considering a multi- 
component dark matter sector. From our study of the 
background, we have shown that multiple coupled flu- 
ids can mimic a single coupled species. In particular, 
multiple fluids with constant couplings (well motivated 
in particle physics) can behave like a single fluid with 
a time varying coupling. Additionally, the cosmological 
background can behave like a ACDM universe during the 
matter dominated epoch, for models in which the effec- 
tive coupling vanishes. 

We established that perturbations allow us to break 
this degeneracy between a single effective dark matter 
species and multiple fluids. Specifically, observers could 
distinguish between a single fluid with a time varying 
coupling /?(</>) and an effective fluid with an apparent 
time varying coupling. 

Perturbation growth is affected on all scales by the 
couplings. On scales smaller than the horizon, particles 
feel a new force, mediated by the scalar field. Depend- 
ing on the relative sign of the coupling, this force can 
be either attractive or repulsive. Depending on the cou- 
pling strength and relative densities, the growth of per- 
turbations for a given species can be either enhanced or 
suppressed. We remark that the instabilities observed in 
[53, 54] are not present in the models considered here. 

Compared to previous studies with a single fluid, CMB 
constraints on the strength of multiple couplings can be 
relaxed. We find that the system is not constrained to 
have Pi < 0.1. We expect stronger constraints to come 
from large scale structure surveys probing the matter 
power spectrum. Additionally, supernovae observations 
will constrain the evolution of the apparent equation of 
state, which can significantly deviate from w = — 1 and 
can even be w < — 1. Observational constraints on these 
models should be considered in future work. Addition- 
ally, a study of non-linear perturbations on small scales 
should be undertaken, to identify further signatures of 
these interactions (such as in [43] and [55]). For example, 
the new interactions could modify the density profiles of 
dark haloes or induce a bias between baryons and dark 
matter. For a single coupled species, this was studied in 
[55] , who conclude that cosmologies with dark-dark cou- 
plings are viable from an observational standpoint. It 
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would be interesting to generalise these findings to our The authors acknowledge support from STFC. 
models, but this is outside the scope of this paper. 
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